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We investigate here the genericity and stability aspects for naked singularities and black holes that arise as 
the final states for a complete gravitational collapse of a spherical massive matter cloud. The form of the matter 
considered is a general Type I matter field, which includes most of the physically reasonable matter fields such 
as dust, perfect fluids and such other physically interesting forms of matter widely used in gravitation theory. We 
first study here in some detail the effects of small pressure perturbations in an otherwise pressure-free collapse 
scenario, and examine how a collapse evolution that was going to the black hole endstate would be modified and 
go to a naked singularity, once small pressures are introduced in the initial data. This allows us to understand 
the distribution of black holes and naked singularities in the initial data space. Collapse is examined in terms of 
the evolutions allowed by Einstein equations, under suitable physical conditions and as evolving from a regular 
initial data. We then show that both black holes and naked singularities are generic outcomes of a complete 
collapse, when genericity is defined in a suitable sense in an appropriate space. 
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The aim of the present work is twofold. Firstly, we study here the effect of introducing small pressure perturbations in an 
otherwise pressure-free gravitational collapse which was to terminate in a black hole final state. For such a purpose, spherically 
symmetric models of black hole and naked singularity formation for a general matter field are considered, which undergo a 
complete gravitational collapse under reasonable physical conditions while satisfying suitable energy conditions. Secondly, we 
investigate the genericity and stability aspects of the occurrence of naked singularities and black holes as collapse endstates. The 

$_i \ analysis of pressure perturbations in known collapse models, inhomogeneous but otherwise pressure-free, shows how collapse 

&0 final states in terms of black hole or naked singularity are affected and altered. This allows us to examine in general how generic 

these outcomes are and we study in the initial data space the set of conditions that lead the collapse to a naked singularity and 

p^ . investigate how 'abundant' these are. While it is known now for some time that both black holes and naked singularities do arise 

£> ' as collapse endstates under reasonable physical conditions, this helps us understand and analyze in a clear manner the genericity 
0^ . aspects of occurrence of these objects in a complete gravitational collapse of a massive matter cloud in general relativity. 

The physics that is accepted today as the backbone of the general mechanism describing the formation of black holes as the 
endstate of collapse relies on the very simple and widely studied Oppenheimer-Snyder-Datt (OSD) dust model, which describes 

• ' ■ the collapse of a spherical cloud of homogeneous dust \m, yfl. In the OSD case, all matter falls into the singularity at the same 
comoving time while an horizon forms earlier than the singularity, thus covering it. A black hole results as the endstate of 
collapse. Still, homogeneous dust is a highly idealized and unphysical model of matter. Taking into account inhomogeneities 
in the initial density profile it is possible to show that the behaviour of the horizon can change drastically, thus leaving two 
different outcomes as the possible result of generic dust collapse: the black hole, in which the horizon forms at a time anteceding 
the singularity, and the naked singularity, in which the horizon is delayed thus allowing null geodesies to escape the central 
singularity where the density and curvatures diverge, to reach faraway observers [3]-|7]. 

It is known now that naked singularities do arise as a general feature in General Relativity under a wide variety of circum- 
stances. Many examples of singular spacetimes can be found, but their relevance in models describing physically viable scenarios 
has been a matter of much debate since the first formulation of the Cosmic Censorship Hypothesis (CCH) [8]. In particular, the 
formation of naked singularities in dynamical collapse solutions of Einstein field equations remains a much discussed problem 
of contemporary relativity. The CCH, which states that any singularity occurring in the universe must be hidden within an event 
horizon and therefore not visible to faraway observers, has remained at the stage of a conjecture for more than four decades now. 
This is also because of the difficulties lying in a concrete and definitive formulation of the conjecture itself. While no proof or 
any mathematically rigorous formulation of the same exists in the context of dynamical gravitational collapse (while some proofs 
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exist for particular classes of spacetimes that do not describe gravitational collapse, as in [9] and ifToil ). many counterexamples 
have been found over the past couple of decades flllll - irifll . Many of these collapse scenarios are restricted by some simplifying 
assumptions such as the absence of pressures (dust models) or the presence of only tangential pressures lfl7ll - ll23ll . It is well- 
known that the pressures cannot be neglected in realistic models describing stars in equilibrium. It seems natural therefore that 
if one wishes to study analytically what happens during the last stages of the life of a massive star, when its core collapses under 
its own gravity thus forming a compact object as a remnant, pressures must be taken into account. 

Therefore, further to early works that showed the occurrence of naked singularities in dust collapse, much effort has been 
devoted to understanding the role played by pressures [24]-[30]. The presence of pressures is a crucial element towards the 
description of realistic sources as we know that stars and compact objects are generally sustained by matter with strong stresses 
(either isotropic or anisotropic). At first it was believed that the naked singularity scenario could be removed by the introduction 
of pressures, thus implying that more realistic matter models would lead only to the formation of a black hole. We now know 
that this is not the case. The final outcome of collapse with pressure is entirely decided by its initial configuration and allowed 
dynamical evolutions and it can be either a black hole or a naked singularity. Furthermore it is now clear that within spherical 
collapse models (be it dust, tangential pressure or others) the data set leading to naked singularities is not a subset of 'zero 
measure' of the set of all possible initial data. 

Despite all this work we can still say that much more is to be understood about the role that general pressures play during 
the final stages of collapse. Perfect fluid collapse has been studied mostly under some simplifying assumptions and restrictions 
in order to gain an understanding, but a general formalism for perfect fluids described by a physically valid equation of state is 
still lacking due to the intrinsic difficulties arising from Einstein equations. Considering both radial and tangential pressures is a 
fundamental step in order to better understand what happens in the ultra-dense regions that forms at the center of the collapsing 
cloud prior to the formation of the singularity. For this reason, perfect fluids appear as a natural choice since these are the models 
that are commonly used to describe gravitating stars in equilibrium and since it can be shown that near the center of the cloud 
regularity implies that matter must behave like a perfect fluid. 

In the present paper we use a general formalism to analyze the structure of collapse in the presence of perfect fluid pressures. 
This helps to understand better realistic collapse scenarios and their outcomes and brings out clearly the role played by pressures 
towards the formation of black holes or naked singularities as the endstate of collapse. We examine what are the key features 
that determine the final outcome of collapse in terms of a black hole or a naked singularity when perfect fluids, without any 
restriction imposed by the choice of an equation of state, are considered. The reason we do not assume an explicit equation of 
state here is that the behavior of matter in ultra-dense states in the final stages of collapse is unknown. On the other hand, having 
regularity and energy conditions satisfied provides a physically reasonable framework to study collapse endstates. 

We find that not only naked singularities are not ruled out in perfect fluid collapse scenarios but also that the separation 
between the black hole region and the naked singularity region in the space of all possible evolutions has some interesting 
features. In particular, we show that the introduction of small pressures can drastically change the final fate of the well-known 
pressureless models. For example, we see that adding a small pressure perturbation to an inhomogeneous dust model leading to 
a black hole can be enough to change the outcome of collapse to a naked singularity, and viceversa. 

Further to this, we investigate here the space of initial data and collapse evolutions in generality, in order to examine the 
genericity of naked singularities in collapse. To study small pressure perturbations as well as the genericity and stability aspects, 
we use the general formalism for spherically symmetric collapse developed earlier 13 ill . H32I1 in order to address the basic 
problem of how generic is a given collapse scenario which leads to the formation of naked singularities. Given the existence 
of an increasing number of models describing collapse leading to a naked singularity, the issue of genericity and stability of 
such models in the space of initial data has become the crucial ingredient in order to decide whether the Cosmic Censorship 
Hypothesis in its present form should be conserved, modified or dropped altogether. 

It should be noted, however, that the concepts such as genericity and stability are far from well-defined in a unique manner in 
general relativity, as opposed to the Newtonian gravity. A major difficulty towards such a task is the non-uniqueness of topology, 
or the concept of 'nearness' itself in a given spacetime geometry B33ll . One could define topology on a space of spacetime 
metrics by requiring that the metric component values are 'nearby' or also additionally requiring that their n-th derivatives are 
also nearby, and in each case the resulting topologies will be different. This is in fact connected in a way with the basic problem 
in arriving at a well-formulated statement of the cosmic censorship itself. There have been attempts in the past to examine the 
genericity and stability of naked singularities in special cases. For example, in [34] it was shown that for certain classes of 
massless scalar field collapse the initial data leading to naked singularity has, in a certain sense, apositive codimension, and so 
the occurrence of naked singularity is unstable in that sense. On the other hand, it was shown in B35I1 . ([36] that naked singularity 
occurrence is stable in the sense of the data sets leading to the same being open in the space of initial data. But these need not 
be dense in this space and so 'non-generic' if we use the definition of 'genericity' in the sense given in the dynamical systems 
theory (where a set of initial data leading to a certain outcome is said to be generic if it is open and dense within the set of all 
initial data). In that case, however, both black hole and naked singularity final states turn out to be 'non-generic'. 

Therefore, in the following we adopt a more physical definition of 'genericity' in the sense of 'abundance', and we call generic 
an initial data set that has a non-zero measure, and which is open (though not necessarily dense) in the set of all initial data. With 
this definition, the results obtained in IB5ll and IBol would mean that both black hole and naked singularity are generic collapse 



endstates. We note that we do not deal here with the self-similar models, or scalar fields, which is a somewhat special case. 
Therefore, the issue of genericity and stability of naked singularities in collapse remains wide open, for spherically symmetric as 
well as non-spherical models and for different forms of matter fields. Our consideration here treats in this respect a wide variety 
of physically reasonable matter fields for spherically symmetric gravitational collapse. 

In section [II] the general structure for Einstein equations to study spherical collapse is reviewed and we describe how the 
equations can be integrated thus obtaining the equation of motion for the system. In Section[III]we examine the structure of the 
initial data sets of collapse leading to black hole and naked singularity to gain an insight on genericity of such outcomes for some 
special models and effect of introducing small pressure perturbations is investigated. Section [IV] then considers the genericity 
aspects of the outcomes of collapse with respect to initial data sets. We prove that the initial data sets leading to black holes and 
naked singularities in the space of all initial data sets for perfect fluid collapse are both generic. Section[V]is devoted to a brief 
discussion on equations of state. Finally, in section[Vl]we outline the key features of the above approach and its advantages, and 
point to possible future uses of the same for astrophysical and numerical applications. 

II. DYNAMICAL EVOLUTION OF COLLAPSE 

In this section we summarize and review the key features on spherical gravitational collapse analysis, and also reformulate 
some of the key quantities and equations, especially those relating to the nature and behaviour of the final singularity curve. 
This will be useful in a later section in analyzing the small pressure perturbations in a given collapse scenario, and subsequently 
towards a general analysis of the genericity aspects of the occurrence of naked singularities and black holes as collapse final 
states. 

The regularity conditions and energy conditions that give physically reasonable models are discussed here. The final stages 
of collapse are discussed, evaluating key elements that determine when the outcome will be a black hole or a naked singularity. 
We shall find a function, related to the tangent of outgoing geodesies at the singularity whose sign solely determines the time 
of formation of trapped surfaces in relation with the time of formation of the singularity. We also analyze here the occurrence 
of trapped surfaces during collapse and the possibility that radial null geodesies do escape thus making it visible. We see how 
both features are related to the sign of the above mentioned function, thus obtaining a necessary and sufficient condition for the 
visibility of the singularity. 

A. Einstein equations 

The most general spacetime describing a spherically symmetric collapsing cloud in comoving coordinates r and t depends 
upon three functions u(r, t), ip(r, t) and R(r, t), and takes the form, 

ds 2 = -e 2 ^dt 2 + e 2 ^dr 2 + R(t, r) W . (II.l) 

The energy momentum tensor reads, 

T\ = -p; T; = Pr ; T e g = Tf = pe , (II.2) 

where p is the energy density and p r and pg are the radial and tangential stresses. The metric functions v, ip and R are related to 
the energy-momentum tensor via the Einstein equations that can be written in the form: 

» - -h ■ ,IL3) 

F' 
P = ~^r- , (II.4) 

H R 2 R' ' V 

■■' - oW-ft-tf * (II.5) 



p + Pr R p + p r 

2k = R '^+ R Jj> ( IL6 ) 

F = RQ.-G + H), (II.7) 

where F is the Misner-Sharp mass of the system (representing the amount of matter enclosed in the comoving shell labeled by 
r at the time t) and for convenience we have defined the functions H and G as, 

H = e' 2v{r ' v ^R 2 , G = e- 2 ^ r ' v) R' 2 . (II.8) 



The collapse scenario is obtained by requiring R < and the central 'shell-focusing' singularity is achieved for R = 0, where 
the density and spacetime curvatures blow up. Divergence of p is obtained also whenever R' — 0, thus indicating the presence 
of a 'shell-crossing' singularity. Such singularities are generally believed to be gravitationally weak and do not correspond to 
divergence of curvature scalars, therefore indicating that they are removable by a suitable change of coordinates [37], [38]. For 
this reason in the forthcoming discussion we will be concerned only with the shell-focusing singularity thus assuming R' > 0. 

Since there is a scale invariance degree of freedom we can choose the initial time in such a way so that R(r, ti) = r. Therefore 
we can introduce the scaling function v(r, t) defined by, 

R = rv , (11.9) 

with v(r, ti) = 1, so that the collapse will be described by v < and the singularity will be reached at v = 0. We see that this 
is a better definition for the singularity since at v ^ the energy density does not diverge anywhere on the spacelike surfaces, 
including at the center r = 0. This is seen immediately from the regular behaviour of the mass function near the center that 
imposes that F must go like r 3 close to r — (as it will be shown later). Such a regularity was not clear from equation (1H.4I >. 
especially along the central curve r = 0, where we have R = 0, without the introduction of v. In this manner the divergence 
of p is only reached at the singularity. We note that v acts like a label for successive events near the singularity since it is 
monotonically decreasing in t and therefore can be used as a 'time' coordinate in the place of t itself near the singularity. 

We shall consider the Misner-Sharp mass to be in general a function of the comoving radius r and the comoving time t, 
expressed via the 'temporal' label v as 

F = F(r,v(r,t)) . (11.10) 

Near the center of the cloud this is just equivalent to a rewriting F(r, t). 

It can be shown that vanishing of the pressure gradients near the center of the cloud imply that the radial and tangential stresses 
must assume the same value in the limit of approach to the center [31]. This requirement comes from the fact that the metric 
functions should be at least C 2 at the center of the cloud and is a consequence of the fact that the Einstein equation ( III. 5b for a 
general fluid contains a term in p r — pg that therefore must vanish at r = 0. Since we are interested in the final stages of collapse 
of the core of a star it is therefore straightforward to take that the cloud behaves like a perfect fluid in proximity of r = 0. We 
shall then take, 

Pr =pe =P ■ 

In such a case, we are then left with five equations in the six unknowns p, p, F, v, ip and v. The system becomes closed once 
an equation of state for the fluid matter is defined or assumed, but in general it is possible to study physically valid dynamics 
(namely those satisfying regularity and energy conditions) without assuming a priori an equation of state. In fact it is reasonable 
to suppose that any equation of state that holds at the departure from equilibrium, when the gravitational collapse commences, 
will not continue to hold in the extreme regimes achieved when approaching the singularity. In this case, we are then left with 
the freedom to choose one of the functions. If we take F as the free function then from Einstein equations ( 1II.3I I and ( III .41 . p and 
p will follow immediately and they can be evaluated explicitly once we know v and its derivatives. 

Further, from the requirement for perfect fluid collapse we can write equation ( 111.51 as 

i/=-jL. (hid 

p+p 

Equation ( III. 6b can be integrated once we define a suitable function A(r, v) from 

,J_ 
R 



A tV = v'^-. (11.12) 



Then we get, 

G(r,v) = b(r)e 2A{r - v) , (11.13) 

where the integration function b(r) can be interpreted in analogy with the dust models and is seen to be related to the velocity of 
the infalling matter shells. 

Finally, from equation dll.7l ) we can write the equation of motion for R in the form of an effective potential as 

R 2 = e 2v (- + G-1 j . (11.14) 

This allows us to study the dynamics of the collapsing system in analogy of the usual phase-space tools of the classical mechanics 
type models. 



Equation (III. 141 ) can be expressed in terms of the scaling factor v as, 



F be 2A - 1 

— + 2 > ( IL15 ) 

where the minus sign has been considered in order to study the collapse. We see that in order to have a solution we must have 
' > 0, we can call this a 'reality condition' that is necessary for the collapse dynamics to occur. Solving the 
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equation (111. 1 51 > solves the set of Einstein equations. 

B. Regularity and energy conditions 

Einstein equations provide the relations between the spacetime geometry and the matter distribution within it, however, they do 
not give any statement about the type of matter that is responsible for the geometry. On a physical ground, not every type of matter 
distribution is allowed, and therefore some restrictions on the possible matter models must be made based on considerations of 
physical reasonableness. This usually comes in the form of energy conditions ensuring the positivity of mass-energy density. 

Further, regularity conditions must be imposed in order for the matter fields to be well-behaved at the initial epoch from 
which the collapse evolves and at the center of the cloud. Firstly, the finiteness of the energy density at all times anteceding the 
singularity and regularity of the Misner-Sharp mass in r = imply that in general we must have, 

F(r, t) = r 3 M(r, v(r, t)) , (11.16) 

where M is a regular function going to a finite value Mo in the limit of approach to the center. If F does not go as r 3 or higher 
power in the limit of approach to the center r = 0, we immediately see from the Einstein equation for p that there would be a 
singularity at the center at the initial epoch, which is not allowed by the regularity conditions as we are interested in collapse 
from regular initial configurations. Also, requiring that the energy density has no cusps at the center is reflected in the condition, 

M'(0,v)=0. (11.17) 

As seen before, the behaviour of the pressure gradients near r — suggests that the tangential and radial pressures become 
equal in limit of approach to the center, thus justifying our assumption of a perfect fluid type of matter. Further since the gradient 
of the pressures must vanish at r = 0, we see that p' ~ r near r = which for the metric function v implies that near the center, 

v(r,t)=r 2 g(r,v(r,t))+g (t) , (11.18) 

where the function go(t) can be absorbed in a redefinition of the time coordinate t. From the above, via equation ( III. 121 ) we can 
write 

A,„ = ^>. (11.19) 

From the analogy with the Lemaitre-Tolman-Bondi (LTB) models we can evaluate the regularity requirements for the velocity 
profile b(r). Since near the center b can be written as, 

b(r) = l + r 2 b (r) , (11.20) 

we now see how to interpret the free function b{r) in relation with the known LTB dust models. In fact the cases with bo constant 
are equivalent to the bound {bo < 0), unbound (bo > 0) and marginally bound (bo = 0) LTB collapse models. Also, the condition 
that there be no shell-crossing singularities may imply some further restrictions on b. 

As is known, matter models describing physically realistic sources must be constrained by some energy conditions. The weak 
energy condition in our case implies 

/Q>0, p + p>0 . (11.21) 

The first one is achieved whenever F' > 0. In fact from equation (111.41 ) we see that positivity of p is compatible with positivity 
of R' only if F 1 > 0. Therefore we must have 

3M + rM' > , (11.22) 

that close to the center will be satisfied whenever M(0, v) > 0. Now from M' = M r + M yV v' we can rewrite p as 



3M + rAf r - M <v v 
v 2 R' 



P= 13§ --*■ ( n - 23 ) 



from which we see that the second weak energy condition is satisfied whenever 

3M + rM >r - M. v v > . (11.24) 

Finally the choice of a mass profile satisfying the above equation allows us to rewrite the condition ( M.221 i as 

3M + rM. r - M tV v > -M jV R' , (11.25) 

which is obviously satisfied if the pressure is positive (since it implies M v < 0) but can be satisfied also by some choice of 
negative pressure profiles. 

The dynamical evolution of collapse is entirely determined once the initial data set is given Il39ll - ll42ll . Specifying the initial 
conditions consists in prescribing the values of the three metric functions and of the density and pressure profiles as functions of 
r on the initial time surface given by t = U. This reduces to defining the following functions: 

p{r,U) = pi(r), p(r,U)=pi(r), R(r,U) = Ri(r), u(r,v(r,ti)) = Vi(r), ip(r,ti) =ipi(r). 

At the initial time the choice of the scale function v is such that R L = r, furthermore, from R[ — 1 we get v'(r, t. L ) = 0. 

Since the initial data must obey Einstein equations it follows that not all of the initial value functions can be chosen arbitrarily. 
In fact the choice of the mass profile together with Einstein equations is enough to specify the four remaining functions. From 
equations dH.4b and ( III. 16b . writing 



we get 

while from equation (III.3b we get, 
From equation ( M.51 l we can write, 
with gi (r) related to Mj by 



M^r) = M(r, v(r, U)) = M(r, 1) , (11.26) 

Pi = 3Mi(r) + rMl{r) , (11.27) 

V l = ~{M, v )i . (11.28) 

Vi(r) = r 2 g(r, v(r, U)) = r 2 9l {r) . (11.29) 



2rg l + r 2 g' l = §— . (11.30) 

Pi + Pi 



In turn, V{ can be related to the function A, defined by equation ( III. 12b . at the initial time via equation dll.191 ). 

A, v (r,v(r,U)) = (A, v )i = 2 9l r 2 + g^r 3 . (11.31) 

Finally, the initial condition for if; can be related to the initial value of the function A(r, t) from equation ( III. 6b and equation 

firm . 

A(r,v(r,U)) = Mr) = -^ - \^b(r) ■ 01-32) 

Since we are studying the final stages of collapse and the formation of black holes and naked singularities, we must require 
the initial configuration to be not trapped. This will allow for the formation of trapped surfaces during collapse and therefore we 
must require 



F <(r) _ Ji 



Ri 



r 2 M z {r) < 1 , (11.33) 



from which we see how the choice of the initial matter configuration Mj is related to the initial boundary of the collapsing cloud. 
Some restrictions on the choices of the radial boundary must be made in order not to have trapped surfaces at the initial time. 
This condition is reflected on the initial configuration for G and H since 1 — JJ — G — H and this condition also gives some 
constraints on the initial velocity. In fact to avoid trapped surfaces at the initial time the velocity of the infalling shells must 
satisfy 

\R\ > Vbe A+ » . (11.34) 



We see that the initial velocity of the cloud must always be positive and that the case of equilibrium configuration where R = 
can be taken only in the limit. 

The consideration of a perfect fluid matter model implies that the Misner-Sharp mass is in general not conserved during 
collapse. Therefore the matching with an exterior spherically symmetric solution leads to consider the generalized Vaidya 
spacetime. It can be proven that matching to a generalized Vaidya exterior is always possible when the collapsing cloud is 
taken to have compact support within the boundary taken at r = r&, and the pressure of the matter is assumed to vanish at the 
boundary J43ll - ll46ll . Matching conditions imply continuity of the metric and its first derivatives across the boundary surface. 
Such a matching is in principle always possible but it should be noted that matching conditions together with regularity and 
energy conditions, might impose some further restrictions on the allowed initial configurations. 

In the following we are interested in the local visibility of the central singularity occurring at the end of the collapse. Therefore 
we will restrict our attention to a neighborhood of the central line r — 0. In this case it is easy to see that there always exist a 
neighborhood for which no shell crossing singularities occur. This is seen by the fact that B! = v + rv 1 and therefore, since v > 
and shell crossing singularities are defined by R 1 = 0, we can always fulfill R' > in the vicinity of the center. Furthermore, 
as it was mentioned before, matter behaves like a perfect fluid close to the center. This can be seen from the fact that regularity 
of the metric functions at the initial time requires that v' does not blow up at the regular center. This in turn implies that from 
equation ( III. 5b we must have p r (0, U) — pe(0, U) = 0, and the condition holds for any time t before the singularity. 

C. Collapse final states 

We study now the possible outcomes of collapse evolution. It is known that in general the final fate of the complete collapse of 
the matter cloud will be either a black hole or a naked singularity, depending on the choice of the initial data and the dynamical 
evolutions as allowed by the Einstein equations. 

In order to analyze the final outcome of collapse we shall perform a change of coordinates from (r, t) to (r, v), thus considering 
t = t(r, v). As mentioned earlier this is always possible near the center of the cloud due to the monotonic behaviour of v. In 
this case the derivative of v with respect to r in the (r, t) coordinates shall be considered as a function of the new coordinates, 
v' = w(r, v). We see that regularity at the center of the cloud implies w — > as r —¥ 0. 

We can then consider the metric function v = v{r, v) which is given by equation ( III. 5b . which now becomes, 

, M tVr v + (M vv v - 2M v )w , ,., 

V = ; r it , (11. Jj) 

(3M + rM, r - M tV v)v 
where now R' = v + rw and , r denotes derivative with respect to r in the (r, v) coordinates. This implies, 

Mr v) f 1 M vr v + (M, vv v-2M, v ) W 

A{r ' V) ~J v (3M + rM, r -M, v v)v rdV ' (IL36) 

For the sake of clarity, we may assume here that near the center the mass function M(r, v) can be written as a series as, 

M(r, v) = M (v) + M x {v)r + M 2 {v)r 2 + o(r 4 ) . (11.37) 

As a regularity condition, we take M\ = and A ~ r 2 . The function A can then be written as an expansion and it takes the 
form, 



with the first terms Ai(v) given by 



A{r, v) = A 2 (v)r 2 + A 3 (v)r 6 + o(r 4 ) , (11.38) 



,1 2M 2 , V + (M 0<vv - ^%*) w, r 

M = \ ^ M " dV > (IL39) 

J v 3M - M , v v 

1 pi 6M 3<V + (M ,vv ~ ^"IT 2 " ) W,rr 

A * = 9 / TM M " dV ■ (IL40) 

2 J v 3M - M , v v 

If we restrict our analysis to constant v surfaces then in equation ( III. 40b we can put w and its derivatives to be zero. On the other 
hand if we approach the singularity along a generic curve we cannot neglect the terms in w. 
The equation of motion (III. 15b takes the form 



M be 2A -l _ A „ 

— + 2 , (H.41) 



which can be inverted to give the function t(r, v) that represents the time at which the comoving shell labeled r reaches the event 



' ! , 



t(r, v) = t t + 



M , be 2A -l 

1 5 



zdv . 



(11.42) 



Then the time at which the shell labeled by r becomes singular can be written as a singularity curve as 



t a (r)=t(r,0)=U + 



M , be 2A -l 



zdv 



Regularity ensures that, in general, t(r, v) is at least C 2 near the singularity and therefore can be expanded as, 

t(r, v) = i(0, v) + xi (v)r + X2{v)r 2 + o{r 3 ) , 
with 

1 



t (0, v) = U + 



zdv 



(11.43) 



(11.44) 



(11.45) 



^L + b (0) + 2A 2 



andxi = fH r=0 andx2 = |0lr=o- 

Of course the situations with discontinuities can be analyzed as well with minor technical modifications to the above formal- 
ism. In our case, assuming that t(r, v) can be expanded implies that the first two terms in the expansion of A(r, v) must vanish. 
As seen before vanishing of first term is consistent with the regularity condition for v that follows from the pressure gradients at 
the center, while vanishing of second term implies that M\ must be a constant, which gives, in accordance with the requirement 
that the energy density has no cusps at the center, that M\ = 0. The singularity curve then takes the form 



t s (r) =t + rxi(0) + r 2 X2 {0) + o(r 3 ) , 



(11.46) 



where to = t(0, 0) is the time at which the central shell becomes singular. 

By a simple calculation, retaining for the sake of completeness all the terms in the expansions of M and expanding v as 

v = g2(v)r 2 + ..., we obtain 



Xi(v) 



1 



^ + h 



"01 



2A, 



-dv 



■ A/„ 



(11.47) 



'^oo 



2A 7 



and 



X2{v) 



where we have defined boj = b '(0) 



3(^ + b 01 + 2A 3 ) 2 



<J2 



1^ + 2A 4 + 2A% + 2b 00 A 2 + b 



00^2 T WQ2 



^ + 6 nn + 2A, 



(M^ + b 



'on 



2Ao 



dv 



(11.48) 



D. Trapped surfaces and outgoing null geodesies 

The apparent horizon is the boundary of trapped surfaces which in general is given by, 

gijRiRj - . (11.49) 

In the case of spherical collapse the above equation reduces to G — H = 0, which together with equation dll.71 gives, 

1 - | = . (11.50) 

This describes a curve v a h{r) given by 

v ah = r 2 M{r,v ah ). (11.51) 
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Inversely, the apparent horizon curve can be expressed as the curve t a h (r) which gives the time at which the shell labeled by r 
becomes trapped. 

In the dust case, the condition M = M(r) implies that approaching the singularity the radius of the apparent horizon must 
shrink to zero thus leaving to as the only point of the singularity curve that can in principle be visible to far away observers. On 
the other hand, in the perfect fluid case we note that models where the mass profile has different dependence on v will lead to 
totally different structures for the apparent horizon and the trapped region. Indeed in full generality there can be cases where 
non-central singularities become visible. This is possible in the case in which M{r, v) goes to zero as v goes to zero, leaving — 
bounded (see e.g. [47]). 

Presently, we are interested in the case where only the central singularity would be visible. In order to understand what are 
the features relevant towards determining the visibility of the singularity to external observers we can evaluate the time curve of 
the apparent horizon in such cases as 

t ah (r) = t s (r) - / -.dv , (11.52) 

M , be 2A -l 
v r 2 

where t s (r) is the singularity time curve, whose initial point is to = t s (0). Near r = 0, equation ( lH.52t can be written in the 
form, 

t ah (r) =t + Xi{Vah)r + X2(v ah )r 2 + o(r 3 ) , (11.53) 

from which we see how the presence of pressures affects the time of the formation of the apparent horizon. In fact, all the initial 
configurations that cause \\ ( or Xi m case mat Xi vanishes) to be positive will cause the apparent horizon curve to be increasing, 
and trapped surfaces to form at a later stage than the singularity, thus leaving the possibility that null geodesies escape from the 
central singularity. By studying the equation for outgoing radial null geodesies it is possible to determine that whenever the 
apparent horizon is increasing in time at the singularity there will be families of outgoing future directed null geodesies that 
reach outside observers from the central singularity, at least locally. 

It can be shown that positivity of the first non-null coefficient x«(0) is a necessary and sufficient condition for the visibility 
of the central singularity 03111 . Nevertheless the scenario of collapse of a cloud composed of perfect fluid offers some more 
intriguing possibilities. In fact we can see from equation (III. 5 Oi l that whenever the mass function F goes to zero as collapse 
evolves it is possible to delay the formation of trapped surfaces in such a way that a portion of the singularity curve t s (r) 
becomes timelike. This in turn leads to the possibility that non-central shells are visible when they become singular [48], thus 
introducing a new scenario that is not possible for dust collapse. It is easy to verify that in order for the mass function to be 
radiated away during the evolution the pressure of the fluid must be negative at some point before the formation of the singularity. 
Despite this seemingly artificial feature negative pressures are worth investigating as they could point to a breakdown of classical 
gravity and could describe the occurrence of quantum effects close to the formation of the singularity [49]. 

III. SMALL PRESSURE PERTURBATIONS OF COLLAPSE 

We will now use the general formalism developed above to study how the outcomes of gravitational collapse, either in terms 
of a black hole or naked singularity, are altered once an arbitrarily small pressure perturbation in the initial data set is introduced. 

The Lemaitre-Tolman-Bondi (LTB) model ( 15011 -15211') for inhomogeneous dust and homogeneous perfect fluid is reviewed 
describing necessary conditions for the LTB collapse scenario. Then certain perfect fluid models are given, using the treatment 
above, by making specific choices for the free function so that it reduces to the LTB case for some values of the parameters. We 
show how the choice of these parameters or introduction of small pressure perturbations is enough to change the final outcome 
of collapse of the inhomogeneous dust. 

From equations (III. 47b we see that if we account for regularity and physically valid density and pressure profiles (typically 
including only quadratic terms in r) we have xi — 0. Then the final outcome of collapse will be decided by xi as written in 
equation tH,48l . We see immediately that once the matter model is fixed globally, thus specifying M, the sign of X2 depends 
continuously on the values at the initial time taken by the parameters Mi, &02 and A, (with j = 2, 4). By continuity then we 
can say that, away from the critical surface for which X2 = 0, if a certain initial configuration leads to a black hole (thus having 
X2 < 0), then changing slightly the values of the initial parameters M, bo, A will not change the final outcome. The same result 
holds for naked singularities and leads us to conclude that every initial data set for which X2 > will have a small neighborhood 
leading to the same outcome [53]. The same, however, cannot be said for the surface separating these two possible outcomes of 
collapse, where X2 = 0. In this case it is the sign of the next non-vanishing x% that determines the final outcome and it is easy to 
see that the introduction of a small pressure such that Xj becomes non-zero for some j < i can change the final fate from black 
hole to naked singularity and viceversa. 

Consider the scenario where the coefficients Xi vanish for every i. This critical surface represents the case of simultaneous 
collapse, or when t s (r) — const, where a black hole forms at the end of collapse and it includes (though it is not uniquely 
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restricted to) the Oppenheimer-Snyder-Datt homogeneous dust collapse model. While this is the case for homogeneous dust, it 
is also easy to show that for inhomogeneous dust and for perfect fluid collapse also there are initial configurations that lead to 
simultaneous collapse, once inhomogeneities, velocity profile, and pressure are chosen suitably. 

In fact if we consider collapse of general type I matter fields, we can always suitably tune the parameters in order to have 
simultaneous collapse and therefore a black hole final outcome. Nevertheless, in all these cases, the introduction of a small 
pressure can drastically change the final outcome by making some Xi turn positive. Of course in full generality there will also be 
regions in which the 'reality condition' is not satisfied and therefore no final outcome is possible. But if we restrict our attention 
to a close neighborhood of the center we will always have a complete collapse of the inner shells thus leading to a black hole or 
a naked singularity. 

In this sense we can consider a small perturbation of any type I fluid collapse. We see that the initial data not lying on the 
critical surface will not change the outcome of collapse once a small perturbation in M or p or b is introduced. On the other 
hand, those initial data sets that belong to the critical surface might indeed change outcome entirely as a result of the introduction 
of a small inhomogeneity, or a small pressure or small velocity. We shall now consider below some collapse models that can be 
obtained from the above formalism, and analyze these under the introduction of small pressure perturbations. 

A. Lemaitre-Tolman-Bondi collapse 

The simplest model that can be studied for small pressure perturbations is the well-known Lemaitre-Tolman-Bondi spacetime, 
where the matter form is dust with pressures assumed to be vanishing. It is interesting to know how the collapse outcome 
would change when small pressure perturbations are introduced in the cloud, which is a more realistic scenario compared to 
pressureless dust. The spacetime metric in this case takes the form, 



ds 2 = -dt 2 



Rf 



1 



r2 M 7 



-dr 2 + R 2 dfl 2 



(III.l) 



and it can be obtained from the above formalism not only if we impose the matter to be dust but also once we require homogeneity 
of the pressures (namely imposing p' = 0). In fact, if we take p = p(t) or p — 0, from Einstein equations, together with the 

regularity condition for v and b(r) we obtains = OandG = l+r 2 bo(r). The equation of motion (III. 151 ) becomes v — y-^f + &o 
where M is a function of t only in the case of homogeneous pressure, and it is a function of r only in the case of dust. In the case 
where M = Mq is a constant we retrieve the Oppenheimer-Snyder-Datt homogeneous dust model that, as it is known, develops 
a black hole at the end of collapse. 

The inhomogeneous dust model is obtained by requiring M — M(r). In this case from equation ( III. 31 ) follows p = and in 
general v can be a function of r and t (requiring v = v{t) is a necessary and sufficient condition to obtain the OSD case). The 
final outcome of collapse is fully determined once the mass profile and the velocity profile are assigned (see figure[T]i. 
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0,4- 




BH 


0,2- 






o- 




\ (0,0)= OSD 


0,2- 


NS 




0,4- 







-0,2 



-0,1 



M, 



FIG. 1: The phase space of initial data for the LTB collapse model with Mo = 1 and M\ — 6oo = &oi = 0. In this case \i = and \2 
determines the final outcome of collapse depending on the values of M2 and 602 ■ There are initial data sets that have a whole neighborhood 
leading to the same outcome in terms of either a black hole or naked singularity. The OSD case lies on the critical surface separating the two 
outcomes which is defined by Xi = for all i. 



11 

For example, in the marginally bound case (namely bo = 0) the singularity curve for inhomogeneous dust becomes t s (r) — 
ti + 1—7^ — and the apparent horizon curve is given by t a h{r) = t s (r) — |r 2 M(r) and in general collapse may lead to 

y M(r) 

black hole or naked singularity depending on the behaviour of the mass profile M(r). By writing M(r) near r = as a series 
M(r) = Mq + M\r + M 2 r 2 + o(r 3 ) we see that the lowest order non-vanishing M, (with i > 0) governs the final outcome of 
collapse. As expected in this case, we have 

Xi (0) = -l^r, (HI.2) 



3 M* 
1 M\ 


1 M 2 


4 M? 


3 M^ 



X2(0) = j^-^—t- (HO) 

4 M 2 d M 2 

In the perfect fluid LTB model (corresponding to the FRW cosmological models in case of expansion) we have that requiring 
v = v(t) is a sufficient condition for having homogeneous collapse. In fact the following two statement can be easily proved: 



1. v = v(t) & 



2. v = v(t) 



\M(r,t)=M(t) 
[bo(r) = k . 

\p = p{t), 
\p = p(t). 



The overall behaviour of the collapsing cloud is determined by the three functions M(r, t), v(r, t) and &o( r ) (which, as we 
have seen, are not independent from one another) and the special cases of Oppenheimer-Snyder-Datt metric and Lemaitre- 
Tolman-Bondi perfect fluid metric can be summarized as follows. 

Mo, k, v(t) <^> Homogeneous dust 

M(r), bo(r), v(r,t) <=> Inhomogeneous dust 
M(t), k, v(t) => Homogeneous perfect fluid 

B. Perturbation of inhomogeneous dust 

We consider now an example based on the above framework by introducing a small pressure perturbation to the inhomoge- 
neous dust model described in the previous section. We consider in general v = v(r, t) and the mass function is chosen of the 
form 

M = Mo + M 2 (v)r 2 , (III.4) 

where Mq is a constant and the pressure perturbation is taken to be small in the sense that Mq 3> \M$ \ at all times (in this way, 
as collapse progresses and the density diverges the model remains close to the LTB collapse as the pressure is always smaller 
than the density). We immediately see that setting M 2 — C reduces the model to that of inhomogeneous dust (and C = further 
gives the OSD homogeneous dust). We note that in this case no non-central singularities are visible since M does not vanish at 
v = 0. Therefore, just like in the dust case, only the central singularity at r = might eventually be visible. 
In this case we can integrate equation ( |H.39t explicitly to obtain 

where C now is the value taken by M 2 at the initial time when v = 1. Therefore we can take the mass function in such a way 
that it corresponds to the LTB inhomogeneous dust at initial time with the pressure perturbation being triggered only at a later 
stage. We can therefore take, 

M 2 (v) = C + e(v) , (111.6) 

and the initial condition M<z(l) = C implies that e(l) = (remember that v € [0, 1]). By taking all the higher order terms to be 
vanishing we easily see that Ai = for i > 2 and 



12 



near r = 0. 

In this case the pressure and the energy density near r = become 



P 



' 2 

-r , p 



3M 



.C- 



(III.8) 



We therefore have two simple possibilities for the choice of the free function e which determines M: 

1. e > which implies e :V < and positive pressure. 

2. e < which implies e „ > and negative pressure. 

From the above, further assuming bo = for simplicity in accordance with marginally bound LTB models, it follows imme- 
diately that xi(0) = and 



X2(0) = -- 



1 f 1 (C + H(v))y/v~ 



(m 



dv 



4 ev 
3 M 



9M 2 j 



1 ev(e + e jV v)s/v 



dv 



(111.9) 



M 



4 ev 
3 Mo 



where we defined the function 



H(v) 



r , v v 



(III. 10) 



We see that \2 is divided into two integrals. If we assume that the pressure perturbation is small (thus considering Mo to be 
big) then the second integral can in principle be neglected. In fact for a suitable choice of Mo it is not difficult to prove that 
the function at the denominator will be positive and monotonically increasing, and therefore it would not affect the sign of the 
integral, while the second integral will be small enough as compared to the first one. Positivity of \i w iU men t> e decided by the 

signofC + ff(w). 
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FIG. 2: The phase space of initial data for perfect fluid collapse with il/o = 1, bo — and M = Mo + (C + e(v))r' 2 . On the left 
the pressure is taken to be p = e(3 — - — -^)r 2 and it reduces to zero for s = 0. On the right the pressure perturbation is given by 
p{r,v) = (0.3 + 2j2j/— 2^-)r 2 . The introduction of the pressure can uncover an otherwise clothed singularity depending on the values of e, 
y and C. Different choices of M with e(v) = t/o + J/i« + V2V 2 + y3V 3 will have a similar qualitative behaviour. 



In order to have naked singularity we must have \2 > 0. This is certainly the case whenever 

C < -H(v) , 



(III. 11) 



for any v. Therefore if we define D\ = min{iJ(u), v <G [0, 1]} all the values of C < D\ will lead to a naked singularity. On 
the other hand, values of C > D^ = ma.x{H(v), v € [0, 1]} will lead to the formation of a black hole. For C £ [D\, D2} the 
explicit form of e(v) is what determines the sign of \2- It can be proven that we can have models in which positive values of C 
lead to the formation of a naked singularity whereas the inhomogeneous dust case led to a black hole. 
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C. A perturbation with pressure p(r) 

We now analyze another case where the pressure perturbation introduced does not depend explicitly on v. A similar situation 
was investigated by one of us earlier in B35I1 . Here we consider a pressure perturbation of LTB models of the form p = p(r). 
If we impose 

M tVV v - 2M, V = , (111.12) 

in the equation dH.351 >. we can then solve for M and explicitly evaluate G. We obtain in that case, 

M(r, v) = y(r)v 3 + z{r) , (111.13) 

where y{r) and z{r) come from integration. Here the case where y — reduces the system to dust. In this sense, if we keep 
\y\ -C z we will consider this model to be a small perturbation of LTB in a similar way as was discussed in the previous example. 
Then the radial derivatives of z will correspond to the inhomogeneities in the LTB models while the pressure, or the function 
y, can be taken to be either positive or negative (with positive pressure corresponding to an increasing mass function M and 
negative pressure corresponding to a decreasing mass function). 

In order to work on a specific model for the sake of clarity, we assume that y and z can be expanded near the center as, 

y(r) = ya + yir + y 2 r 2 + ..., (111.14) 

z(r) = z Q + z ir + z 2 r 2 + ... . (111.15) 



Then the pressure and density become 



q i \ 3z + z, r r + y, r rv 3 

p = -3y(r),p= — p. (III. 16) 

v z R' 



Imposing regularity requires that y, r (0) = z, r (0) = 0, which implies y\ = Z\ = 0. At the center of the cloud the pressure and 
density become po = — 3yo an d Po — 3 (■% + Ho), an d the energy conditions impose that zq > and j/q > ~ z o- Then from 
equation ( III. 36b we can integrate explicitly to obtain, 

A(r,v)=]n( U ^ + \ ) , (III. 17) 

\u{r) + v i J 



where we have defined 



u(r) = 3z(r) + Z - r . (111.18) 



y, r r 

Then we get, 

G(r,v)=b(r)(^±\) . (III. 19) 

From the above expressions, we can easily obtain now \\ an d Xi- F° r simplicity we consider here the case where 6o = ?/4 = 0. 
Then \i = and we get, 



1 
X2 = ~2 



2^v 6 V2V 3 + z 2 + l^v{v 3 - 1) [y 2 (v 3 - |) + f z 2 ] 

/ I 

^/yov* + z + |f w(l - v 3 ) ^ yoV 3 + ZQ + ±Ul v {1 - w 3 )) ' 



y/vdv (111.20) 



We see from here that the sign of \ 2 is explicitly determined by the inhomogeneities (z 2 ) and the pressure gradient (y 2 ) (see 
FigOJ. 

Once again taking yo = y 2 = and z 2 = reduces the system to the OSD collapse scenario and we see that the introduction 
of the slightest pressure can change drastically the outcome of collapse. On the other hand, taking only y$ = y 2 = (with 
z 2 y^ 0) we retrieve the LTB model and once again to change the final outcome of collapse we must choose y 2 suitably to 
balance the contribution to xi given by the inhomogeneities. Therefore we see again that also within this perturbation model any 
collapse with initial data taken in a neighborhood of a model leading to a certain outcome and not lying on the critical surface 
will result in the same endstate. 



The equation for the apparent horizon curve can be easily written in this case and becomes 

r 2 y(r)v ah (r) 3 - v ah (r) + r 2 z(r) = , 
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(111.21) 



which is a cubic equation in v a h that admits in general three solutions in the case where 27r z f- > 0. Obviously this 

condition is satisfied near r — for positive pressures (that correspond to negative y) and this indicates, as already stated, that 
the mass function is not vanishing at any time and therefore the central shell becomes trapped at the time of formation of the 
singularity. On the other hand, for negative pressures the central shell is not trapped and the formation of the apparent horizon 
can be shifted to some outer shells or removed altogether. 




-0,08 -0,06 -0,04 -0,02 



FIG. 3: The numerical evaluation of the sign of X2 as a function of yi and zi provides the phase space of initial data for the perfect fluid model 
with p = p(r). Here we have taken zo = 1 and j/o = —1/2. Again, there are initial data sets that have a whole neighborhood leading to 
the same outcome which is either a black hole or naked singularity. The OSD case lies on the critical surface separating the black holes from 
naked singularities. 



D. Perturbation of a general simultaneous collapse 



As we mentioned, the case of simultaneous collapse, which means that the final state of collapse is necessarily a black hole, 
need not be restricted to the Oppenheimer-Snyder-Datt model only. In fact for different kinds of general type I matter fields 
there might be suitable choices of the parameters that lead the final state of collapse to be simultaneous. This is of course the 
case when the pressures are homogeneous, that is represented by the time reversal of the Friedman-Robertson-Walker model, 
but more general matter models might also lead to the same behaviour. 

Simultaneous collapse means that all matter shells terminate into the singularity at the same time. Then we see from equation 
( 111.46b , which describes the singularity curve, that all coefficients Xi(0) must vanish, or equivalently that t s (r) = to- From 
equation ( 1II.431 I we can see that a sufficient condition for simultaneous collapse is 



-2v 



M 



be 2A -l 



h{v). 



(111.22) 



This condition, for any given choice of the matter model, leads to a choice of the free function b(r) as 



b(r) 



e- 2A [ 1 



-2// 



M 



It is easy to check that in the case of dust this reduces to 



bo(r) 



M(r) 



h(v) 



(111.23) 



(111.24) 
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and since the mass profile in this case is a function of r only, we conclude that equation ( 1III.221 I is satisfied for dust only by 
the case of homogeneous dust collapse (where M(r) = Mq and bo — k). Nevertheless, as we have said, this need not be the 
only case when the collapse is simultaneous. For example it's straightforward to see that when pressures are considered, the 
same condition as above holds for collapse of an homogeneous perfect fluid, where, in this case, M = M(t). Furthermore it is 
possible that the condition ( IIII.22I ) can be satisfied by some suitable function b(r) also for more general pressure profiles, since 
in general the mass function depends on both r and t, or for some other suitable choice of the velocity profile. 

In order to better understand the conditions under which we can have simultaneous collapse in full generality, let us now 
consider a general perfect fluid matter model given by a choice of M(r,v) (which implies A from equation ( III. 361 )). thus 
specifying all coefficients Mi(v). Firstly, we notice from equation ( III. 451 ) that a suitable choice of &oo is necessary in order 
for the reality condition to be fulfilled near the center. Therefore, once we made this choice and carried out the integration for 
equation ( III. 47b we see that \i (0) depends linearly on boi only. In fact we can write 



Xi(0) = ~2 / ai(v)dv--j- 



Pi(v)dv , 



(111.25) 



from which we see that it will always be possible to choose &oi suitably such that xi(0) = 0. The same reasoning can then be 
applied for all other coefficients Xi(0) that will depend linearly on b^i as, 



XiW 



i ,i 

ai(v)dv + b 0i / (3 t (v)dv 



(111.26) 



and therefore for a given mass profile M (r, v) we can have simultaneous collapse if a suitable velocity profile b(r) given by 



6(r)-l+r»5:^ 



(111.27) 



exists. This means that the power series X^o % ir * should converge to some function 6o (r) with a radius of convergence greater 
than the boundary of the cloud. 

This is certainly possible in the case of homogeneous pressures, where the condition that all Xi(0) vanish imposes that 
~l2iLo ^f T% = r5' anc ' therefore b(r) = const. Also, as we have seen in the examples above, this might be possible for other 
matter profiles as well. Given any such model leading to simultaneous collapse (and thus to the formation of a black hole), we 
have shown that the introduction of the slightest pressure perturbation in the initial data can turn the final outcome into a naked 
singularity. 

Overall we have seen that the sign of \i an d therefore the final outcome of collapse shares similar qualitative behaviour in 
different perfect fluid models as it is summarized in Fig. [4] 
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FIG. 4: Illustrative representation of the phase spaces of initial data for different perfect fluid collapse models with assigned mass profiles 
M(r, v) and with bo — 0. The introduction of pressures p and inhomogeneities C can uncover an otherwise clothed singularity. Also, 
different choices of M can have the opposite effect, that is, a singularity that was naked can become covered. The initial data set h has a whole 
neighborhood lying in the initial data space leading to naked singularity and can therefore be considered stable. The same holds for the initial 
data set h leading to black hole. The LTB model (obtained for p = and) can lead to either of the outcomes and the OSD scenario (obtained 
when there are no inhomogeneities) always lies on the critical surface separating the two outcomes and is therefore 'unstable'. 
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IV. GENERICITY OF BLACK HOLES AND NAKED SINGULARITIES 

As we have seen, the final outcome of collapse depends upon the evolution of the pressure, the density and the velocity 
profiles. If the system is closed, as it is in the case where an equation of state describing the relation between p and p is given, 
then specifying the values of the above quantities at the initial time uniquely determines the final outcome of collapse. If the 
system is not closed, then we must further specify the behaviour of the free functions. Once again for every possible choice of 
the free function(s) the final outcome of collapse is decided by the initial values of p, p and b. 

The genericity is defined here as every point in the initial data set leading to a naked singularity (or a black hole) has a 
neighborhood in the space of initial data for collapse whose points all lead to the same outcome. We show that the initial data set 
leading the collapse to a naked singularity forms an open subset of a suitable function space comprising of the initial data, with 
respect to an appropriate norm which makes the function space an infinite dimensional Banach space. The measure theoretic 
aspects of this open set are considered and we argue that a suitable well-defined measure of this set must be strictly positive. 
This ensures genericity of initial data in a well-defined manner. 

At this point, the question of whether the given outcome is 'generic' or not in a certain suitable sense yet to be defined, with 
respect to the allowed initial data sets, arises naturally. We shall therefore analyze the expression ( III. 471 ) for the genericity of 
initial data leading the collapse to a naked singularity. Similar conclusions apply to the case where \i vanishes and we must 
analyze equation ( III. 481 ) and they can be used to investigate the genericity of the black hole formation scenario just as well. 

As is known the concept of 'genericity' is not well-defined in General Relativity. Normally, by the word 'generic', one 
means 'in abundance' or 'substantially big'. This terminology has been used by many researcher working in relativity, and in 
gravitational collapse in particular (see for example, J5411 - ll5aO . In the theory of Dynamical Systems, however, the definition 
of 'genericity' is given more tightly. There, one considers the class V of all C r vector fields (dynamical systems) defined on 
a given manifold. A property P satisfied by a vector field X in V is called generic if the set of all vector fields satisfying this 
property contains an open and dense subset of V. This was the definition used by one of us in IB5I1 and 13611 . However, such a 
definition would render both black holes and naked singularities to be 'non-generic', as we remarked earlier. Therefore in the 
present paper we have opted for a less stringent but physically more meaningful definition of genericity by requiring that the 
subset is open, and that it has a non-zero measure. The main reason for this comes from the fact that the 'denseness' property 
depends on the topology used and the parent space used, and there are no unique and unambiguous definitions available in this 
regard as discussed earlier. Hence the nomenclature of 'generic' in the present paper is used in the sense in which most of the 
relativists use it, i.e. in the sense of abundance. This looks physically more satisfactory definition, allowing both black holes and 
naked singularities to be generic. In any case, the key point is that regardless of the definition used, both the collapse outcomes 
do share the same 'genericity' properties, which is what our work here shows. 

A. Existence of the set of initial data leading to naked singularity 

First of all we note that the functions must satisfy the 'reality condition' for the gravitational collapse to take place, namely 

^l + b {0)+2a(0,v)>0, (IV.l) 

v 

where we have defined a(r, v) = 2 . 

We shall now prove that, given a mass function M(r, v), and the function a(r, v) on the initial surface, there exists a large 
class of velocity distribution functions 60 (r) such that the final outcome is a naked singularity. We choose 60 (r) to satisfy the 
following differential equation on a constant u-surface, 

1 ^f^- + b' (r)+2a'(r,v) , N 

s ou y_^j_ = B{r ^ y) _ (IV2) 

[^+b (r)+2a(r,v) 
for < r < ft, where B(r, v) is a continuous function defined on a domain D = [0, r&] x [0, 1] such that 

M^l jVq) + 2a '( ,v) 

B(0,v) = - 2 3U K -^ < , (IV.3) 



2 



M(0,v) 



^l + b (0) + 2a(0,v) 



for all v e [0, 1]. It will then follow that 

Xi(0) = lim xi{v) = - B(0, v)dv > 
v ->° Jo 



(IV.4) 
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As seen above, this condition ensures that central shell-focusing singularity will be naked. 

We prove the existence of bo(r) as a solution of the differential equation ( IIV.2I I with initial condition ( IIV.3I ) which B(r, v) will 
satisfy. For this purpose, we define 



x(r, v) = M{r, v) + vbo(r) + 2va{r, v) 



Then equation dIV.2b can be written as 



^ [M'(r,v)+vV (r) + 2va'(r,v)] 1 /- f 
2 [M(r,v)+i;&o(r) + 2va(r,v)] 5 2 a; 5 



or 



with the initial condition 



dx 

dr 



2B(r,v)xz , = f(x,r) , 



x{0, v) = M (v) + vb a (0) + 2va(0, v) . 



(IV.5) 



(IV.6) 



(IV.7) 



(IV.8) 



We now ensure the existence of a (^-function x(r, v) as a solution of the above initial value problem defined throughout the 
cloud. The function f(x, r) is continuous in r, with x restricted to a bounded domain. With such domain of r and x, f(x, r) is 
also a (^-function in x which means f(x, r) is Lipschitz continuous in x. Therefore, the differential equation dIV.7b has a unique 
solution satisfying the initial condition ( IIV.8I ). provided f(x, r) satisfies a certain condition given below. 

Further, we can ensure that the solution will be defined over the entire cloud, i.e. for all r € [0, r&], by using the freedom in 
the choice of arbitrary function B(r, v). For this, we consider the domain [0, r&] x [0, d] for some finite d. 

Let us take S = sup |/(a;,r)|. Then the differential equation ( !IV.7t has a unique solution defined over the entire cloud 
provided, 



^ ■ ft d \ d 
n < mf (r 6 , -g) = -g 



(IV.9) 



This condition is to be satisfied according to usual existence theorems to guarantee existence of a unique solution (see for 
example [57]). Equation flIV.9t implies S < — i.e., 



max 

0<r<r b ,0<v<d 



1 



2B(r,v)xi 



< 



>'b 



This, in turn, will be satisfied if 



0< 



B(r,v)x* 



< 



~2^n 



(IV. 10) 



(IV. 11) 



d-Vv . 



for all r € [0, Tf,], 

The collapsing cloud may start with r& small enough so that the expression ^p which is always positive, satisfies the condition 
(1IV.1U with x restricted to a bounded domain. We then have infinitely many choices for the function B(r, v), which is continuous 
and satisfies conditions ( UV31 l and ( lIV.lll i for each choice of v. For each such B(r, v), there will be a unique solution x(r, v) of 
the differential equation (1IV7I >. satisfying initial condition dIV8t . defined over the entire cloud, and in turn, there exists a unique 
function ba(r) for each such choice of B(r, v), that is given by the expression 



&o(r) 



x(r, v) — M(r, v) — 2va{r, v) 



(IV 12) 



Thus, we have shown the following: For a given constant u-surface and given initial data of mass function F(t, r) = r 3 M(r, v) 
and a(r, v) satisfying physically reasonable conditions (expressed on M), there exists infinitely many choices for the function 
b(r) such that the condition (IIV.21 > is satisfied. The condition continues to hold as v -^ 0, because of continuity. Hence, for all 
these configurations the central singularity developed in the collapse is a naked singularity. 

The above analysis shows that the initial data satisfying conditions (II VU . (IIV.3t and (II V lit lead the collapse to a naked 
singularity. If we change the sign in condition (IIV.3t . call it condition ( lIV3l l'. then above analysis apply and the initial data 
satisfying conditions dlV.U . dIV.3l > / and dlV.llt lead the collapse to a black hole. In the cases discussed above, in addition to 
above conditions, energy conditions are also to be satisfied, and we have shown above that this is always possible for matter 
models leading to both possible outcomes. 

Thus, from the above analysis, we get the following conditions which should be satisfied by the initial data in order that the 
end state of collapse is a naked singularity: 
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1. Energy conditions: p > , and p + p > 0. 

2. Reality condition given by equation dl V. 1 i above. 

3. Condition on B(r, v): B(0, v) < for naked singularity and B(0, v) > for black hole. 



4. 0< 



B(r, u)x- 



- 2r b ■ 



For convenience, we denote the function 

C( r , v) = Mh^l + &0 ( r ) + 2a (r, v) . 
v 

Then the reality condition (2) becomes C(0, v) > 0. Assuming this condition, condition (3) will be satisfied if and only if 
> ' + b (0) + 2a' (0, v) < for naked singularity, and > for a black hole. Whenever C(r, v) is an increasing function of 

r in the neighborhood of 0, we get its derivative positive, and so B(0, v) > 0, and end state will be a black hole. On the other 

hand, if C(r, v) is a decreasing function of r in the neighborhood of 0, we get its derivative negative, and so -6(0, v) < 0, and 

the endstate will be a naked singularity. 

Regarding condition (4), using the expression x — vC(r, v), it becomes C'(r, v) < — which will be satisfied if v and Tb are 

sufficiently small. Thus validity of all these conditions does not put any stringent restrictions on the initial data. 

The conclusion then is the following: If the initial data consisting of the mass function M(r, v) and function a(r, v) satisfies 

the above conditions, then there is a large class of velocity functions bo(r) such that end state of collapse is either a black hole 

or a naked singularity, depending on the nature of function C(r, v) as explained above. 

B. Measure of the set of initial data leading to naked singularity 

We now show that the set of initial data Q = {M(r, v), a(r,v), b(r)} satisfying the above conditions which lead the collapse 
to a naked singularity, is an open subset of X x X x X, where X is an infinite dimensional Banach space of all C 1 or C 2 
real-valued functions defined on T> = [0, r&] x [0, d], endowed with the norms 

||Af(r,u)||i = sup|Af|+sup|Af >r |+sup|Af jt) | , (IV. 13) 

V V V 

\\M(r,v)\\ 2 = \\M(r,v) ^ +sup|M. rr | +sup|M. r „| +sup|M.„„| . (IV. 14) 

V V V 

These norms are equivalent to the standard C 1 and C 2 norms, 

|| M(r,v) ||i = sup(|M| + |M,,,| + |M,„|) , (IV 15) 

v 

|| M{r,v) || a = supfll M(r,v) || x | + |M. rr | + \M. rv \ + \M, VV \) . (IV. 16) 

v 

Let Q x = {M(r, v) : M > 0, M is C 1 , E 1 > and E 2 > on V} be a subset of X, where E x = 3M + rM. r + rv'M. v and 
E<z = 3M + rM tr — vM <v , Thus E\ > and E% > are equivalent to energy conditions. 

We first show that Q\ is an open subset of X. For simplicity, we use the C 1 norm, but a similar proof holds for the C 2 norm 
also. For M G Gi, let us put 5 = min(M), 7 = min(£'i) , j3 — min(£ l 2), Ai = max(w), A2 = max(ti') and for r varying in 
[0, Vb] and v € [0, 1], the functions involved herein are all continuous functions defined on a compact domain T> and hence, their 
maxima and minima exist. We define a positive real number 

M= 2 mm{,5 '9'3^'3^'9'3^'3AT } - (IV17) 

Let Mi(r,v) be C 1 in V with || M - M x ||i< p. Using above definition we get \M X - M\ < p, \M lyT - M r | < p, 
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|Mi „ — M„ | < /i over T>. Therefore, for choice of \x, the respective inequalities are 



Mi > M - - > 0, (IV. 18) 



2 



3|Afi-M| < J, 
o 



r\M hr - M, r \ < n\M hr - M r | < |, 
ru'|Afi,„ - M,„| < r 6 A 2 |M M - M,„| < J 



G 



2,\M 1 -M\ < t 
o 

r|Mi, r -M, r | < f, 

u|Mi, w - M, „| < r 6 |Afi, fl - Af,„| < ~, 

D 

that are satisfied on T>. The 2 nd , 3 rd and 4 th inequalities from above yield 

3|Mi -M\ +r\Mi >r -M ir \+rv'\Mi >v -M, v \ < | < 7 < Ej_ . (IV. 19) 

Further, we can write | [3Mi + rMi, r + n/Mi,,,] - £i | < 2<a where £?i > on V. 

Hence, 3Mi + rM\ tT + rv'Mi >v > OonP. Using similar analysis for last four inequalities of equation (1IV. 181 ), we obtain 
3Mi + rMi, r - vMx' v > on V. 

Thus, Mi > 0, Mi is C 1 , [3Mi + rMi, r + rv'Mi t „] > and [3Mi + rMi, r - vM ljV ] > on V provided v > throughout 
T). Therefore, Mi (r, v) also lies in Cq and hence, Q\ is an open subset of X. 

Using the similar argument we can show that the set Q2 = {M(r, v) : \B(r, v)\x^ < -^r-} is also an open subset of X. 
Thus the set of M(r, v) satisfying above conditions forms an open subset of X, since intersection of finite number of open sets 
is open. Similar arguments show that the set of a(r, v) and b(r) satisfying these conditions form separately open subsets of X. 
Hence using definition of product topology we see that the set Q defined above is an open subset of X x X x X. Thus initial data 
leading the collapse to a naked singularity forms an open subset of the Banach space of all possible initial data, and therefore it 
is generic. 

We now discuss measure theoretic properties of the open set Q consisting of the initial data leading the collapse to a naked 
singularity. By referring to the relevant literature about measures on infinite dimensional separable Banach spaces, we argue that 
this Q has strictly positive measure in an appropriate sense. For simplicity, we consider a single space X and its open subset Q. 
We ask the question: Does there exist a measure on X which takes positive value on Ql To answer this question, we note that 
X is an infinite dimensional separable Banach space, and it is a consequence of Riesz lemma in functional analysis that every 
open ball in X contains an infinite disjoint sequence of smaller open balls. So, if we want a translation invariant measure on X 
then its value will be same on each of these balls. Thus, if we demand that the surrounding ball has finite measure, then each of 
these smaller balls will have measure zero. Otherwise sum of their measures would be infinite by countable additivity. In other 
words, for separable Banach spaces, every open set has either measure zero or infinite under a translation invariant measure. So, 
if we wish to have a non-trivial measure on X, then we have to discard the property of translation invariance. In that case we 
must shift to Gaussian measures or Wiener measures. Under these measures, we can conclude that an open subset of X will 
have a positive measure. For example, it is proved in [58] (Theorem 2 on page 159), that Gaussian measure of an open ball in 
a separable Banach space is positive. We can also use Wiener measure on C([0, l])to get the same result (see for example 115911 
andJ60ll). 

However, for all practical purposes in Physics and Astrophysics, physical functions could be assumed to be Taylor expandable. 
Thus assuming that our initial data is regular and Taylor expandable, and again working for simplicity with a single function 
space, instead of product space, we can formulate our problem of measure as follows: Let Y denote the space of Taylor expand- 
able functions defined on an interval [0, T], We consider initial data consisting of functions with first finite number of terms, say 
n terms, which lead the collapse to a naked singularity. These functions will belong to a finite dimensional space isomorphic to 
M. n . Working with supremum norm as above and arguing similarly, we can prove that the initial data set satisfying conditions (1) 
to (4) above is an open subset of R™. Now, we have a standard result (see for example 16111 . prop. 4.3.4, page 83) that a Lebesgue 
measure of an open subset in R n is strictly positive. Denoting this measure by /i„ and the open set by Q n , we get fj, n (Gn) > 0. 
If, further, Q n is bounded, then n n (Q n ) will be finite. Hence normalized Lebesgue measure of an open subset, and in particular, 
of an open ball in R™ is also strictly positive, and in fact bounded. We ask the question: Assuming that fi n {G n ) > 0, can we get 
a measure on Y such that measure of Q is positive? This is answered affirmatively by Maxwell-Poincar theorem which is stated 
as follows (see for example [62]): 
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Consider the sequence of the normalized Lebesgue measures on the Euclidean spheres S 1 ™ -1 C R" of radius r n = Cy/n, c > 
and the limit of spaces 

R 1 CR 2 ... CM" C ... CM 00 . 

Then the weak limit of these measures is the standard Gaussian measure /i which is the infinite product of the identical Gaussian 
measures on the line with zero mean and variance c 2 . Thus the limit exists and is positive on an open subset Q. 

It is also possible to give other approaches which answers affirmatively the existence of such limits which are termed as 
infinite products or in general 'projective limits'. We describe briefly one such approach as described by Yamasaki H63I1 (For 
general concepts on measure theory, we refer to 16411 ). 

Let R°° denote the infinite product of real lines. Let Rg° be the subspace of M°° given by Rg° = {(£n)'- there exists N with 
£„ = for n > N}. Then R°° is the algebraic dual of Rg° and S Rg = is the weak Borel field of R°°. Members of B K g° are called 
weak Borel subsets of R°°. The space Y mentioned above can be seen isomorphic to a subspace of R°°, and is isomorphic to 
Rg° if we consider a finite number of terms in the Taylor expansion. Let m denote one dimensional Lebesgue measure on R. 
Let G = (Gk) be a sequence of Borel sets of R such that < rn(Gk) < oo. We shall define two Borel measures m k and X k by 

m k (B) = m(B)/m(G k ), X k (B) = m k (B n G k ), 

m k is cr-finite, whereas X k is a probability measure on R. 
Consider the product measure 

n oo 

a*g = n mk x n Xk > 

fc=l k=n+l 

then fj,Q ' is cr-finite on R°°. 

Then we have the following theorem: For every weak Borel set E of R°°, put 

n— »oo 

This limit always exists and becomes a a-finite Mg° -invariant measure on E°°. Then \iq lies on 



L(G)=U~ =1 L n (G), 



where 



L n (G)=R n x f[ G k . 

k=n+l 

We note that the measure /j,q defined in this theorem is called the infinite dimensional Lebesgue measure supported by L(G). 
For any open bounded subset E of R°°, E is a Borel set and hence a weak Borel set. Moreover 



H ( G ) (E)=(l[m k x [] X k )(E)=[l[m k )(E)x[ jj ^ ) (E) , 



\ k =\ fc=n+l / \fc=l / \fc=n+l 



and both these factors are finite. Thus the measure [iq takes a non-zero value on an open bounded subset E of R°°. We can 
employ this measure instead of the Gaussian measure mentioned in Maxwell-Poincar theorem to yield the desired result. In any 
case, use of probability measure is inevitable and we conclude that the space of initial data leading to a certain outcome (be it 
black hole or naked singularity), within a specific collapse scenario has non zero measure with respect to the set of all possible 
initial data. 

V. EQUATION OF STATE 

As is known the presence of an equation of state introduces a differential relation for the previously considered free function 
that closes the system of Einstein equations. Examples of simple, astrophysically relevant, linear and polytropic equations of 
states are discussed below. 

In the scenario described above, the relation between the density and pressure could vary during collapse, as it is natural to 
assume in the case where we go from a nearly Newtonian initial state to a final state where a very strong gravitational field is 



21 

present. The equation relating p to p will therefore be represented by some function of r and v that is related to the choice of 
the free function M. There are at present many indications that suggest how in the presence of high gravitational fields gravity 
can act repulsively and pressures can turn negative towards the end of collapse. Therefore if that is the case then the equation of 
state relating pressure and density (which is always positive) must evolve in a non-trivial manner during collapse. 

Typically we can expect an adiabatic behaviour with small adiabatic index at the beginning of collapse when the energy 
density is lower than the nuclear saturation energy. It is not unrealistic to suppose that the equation of state will have sharp 
transitions when matter passes from one regime to another, as is the case when the limit of the nuclear saturation energy is 
exceeded. Towards the end of collapse repulsive forces become relevant thus giving rise to negative pressures and the speed of 
sound approaches the speed of light 116511 . 

Nevertheless it is interesting to analyze the structure of collapse model within one specific regime once a fixed equation of 
state, of astrophysical relevance, is imposed. If we choose the equation of state to be linear barotropic or polytropic we can 
describe collapse of the star right after it departs from the equilibrium configuration where gravity was balanced by the nuclear 
reactions taking places at its center. From an astrophysical point of view, neglecting the energy coming from the nuclear reactions 
occurring at the interior is reasonable since we know that once the nuclear fuel of the star is exhausted the star is subject to its 
own gravity only and the departure from equilibrium occurs in a very short time. In this sense equilibrium models for stars (such 
as the early models studied in the pioneering work by Chandrasekhar [66]) constitute the initial configuration of our collapse 
model and the physical parameters used to construct those equilibrium models will translate into the initial conditions for density 
and pressure. 

As we mentioned before, introducing an equation of state is enough to ensure that the system of Einstein equations is closed 
and so no freedom to specify any function remains. In fact a barotropic equation of state of the form, 

P = Y(p), (V.l) 

introduces a differential equation that must be satisfied by the mass function M(r, v), thus providing the connection between 
equation ( 1H.31 I and equation ( III.4b and making them dependent on R(r, v) and its derivatives only. The dynamics is entirely 
determined by the initial configuration and therefore we see how solving the equation of motion ( III. 1 41 > is enough to solve the 
whole system of equations. 

In this case solving the differential equation for M might prove to be too complicated. Nevertheless with the assumption that 
M can be expanded in a power series as in equation dll.371 ) we can obtain a series of differential equations for each order Mj. 
Expanding the pressure and density near the center we obtain explicitly the differential equations that, if they can be satisfied by 
all Mi converging to a finite mass function M , solve the problem, thus giving the explicit form of M. 

From 

p(r,v) = po{v) + p 2 (v)r 2 + ... = Y (v) + Y 2 (v)r 2 + ... ., 

p(r,v) = po(v)+p 2 (v)r 2 + ... , (V.2) 

with 

Y = Y( Po ) , Y 2 (v) = Y pP2 , (V.3) 

and from Einstein equations ( III. 31 ) and ( III. 41 ) we get 

Y(po) = -% , (V.4) 

v 

Y p p 2 = -% , (V.5) 



with 

3M ° ,vm 

pa = — o- , (V.6) 

5M 2 3M 

Pi = 3 7-W,r{0,V) . (V.7) 

Once again we see that without the knowledge of w, which is related to v', it is impossible to solve the set of differential 
equations in full generality. Furthermore we can see that whenever pressures and density can be expanded in a power series near 
the center the behaviour close to r = approaches that of an homogeneous perfect fluid. 

There are a few equations of state that have been widely studied in equilibrium models for stars and that naturally translate 
into collapse models. The simplest one is a linear equation of state of the form 

P = Ap , (V.8) 
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where A is a constant. This case was studied in B67I1 where it was shown the existence of a solution of the differential equation 
for M, which, from Einstein equations ( III. 3b and ( III. 4b becomes 

3AM + XrM tr + [v + (A + l)n/]Af „ = . (V.9) 

It was shown that both black holes and naked singularities are possible outcomes of collapse depending on the initial data and 
the velocity distribution of the particles. 

Another possibility is given by a polytropic equation of state of the type 

p = \p 1 . (V.10) 

Such an equation of state is often used in models for stars at equilibrium and can describe the relation between p and p in the 
early stages of collapse. Therefore the physical values for po, po, A and 7 at the initial time can be taken from such models at 
equilibrium and expressed in terms of the thermodynamical quantities of the system such as the temperature and the molecular 
weight of the gas. The pressure is typically divided in a matter part (describing an ideal gas) and a radiation part (related to 
the temperature). The exponent 7 is generally written as 7 = 1 + 1/n, where n is called polytropic index of the system and 
is constrained by n < 5 (for n > 5 the cloud has no boundary at equilibrium) 1681 . 16911 . The formalism developed above can 
therefore be used to investigate such realistic scenarios for collapse of massive stars. 

VI. CONCLUDING REMARKS 

We have studied here the general structure of complete gravitational collapse of a sphere composed of perfect fluid without 
a priori requiring an equation of state for the matter constituents, thus allowing for the freedom to choose the mass function 
arbitrarily, as long as physical reasonableness as imposed by regularity and energy conditions is satisfied. 

The interest of such an analysis comes from the fact that the class of perfect fluid models for matter is considered to be 
physically viable for the description of realistic objects in nature such as massive stars and their gravitational collapse. Typically 
perfect fluids are considered to be physically more sound than models where matter is approximated by dust-like behaviour i.e. 
without pressures, or where matter is sustained by only tangential pressures (though the 'Einstein cluster' describing a spherical 
cloud of counter-rotating particles has been shown to have some non-trivial physical validity [70]-[72]). What we have shown is 
that, within the class of perfect fluid collapses, both final outcomes, namely black holes and naked singularities, can be equally 
possible depending on the choice of the initial data and the free function F. In fact our results show that naked singularities and 
black holes are both possible final states of collapse, much in the same way as it has already been proven in the simpler cases 
of inhomogeneous dust and matter exhibiting only tangential stresses. The sets of initial data leading to either of the outcomes 
share the same properties in terms of genericity and stability. 

The structure of initial data sets in the case of OSD, LTB and pressure collapse and their inter-relationship is, however, a 
complicated issue. Nevertheless, we can comment on this based on the studies in 113511 . 13611 . and the results proved in sections Hill 
and|IV]in this paper. In IB611 it was proven that the space of initial data {M (r), b(r)} leading LTB collapse to black hole or naked 
singularity forms an open subset of AT x X, where X is the infinite dimensional Banach space of real C 1 functions defined on the 
domain. As per the analogous results in the case of non-vanishing tangential pressures it follows that the initial data set leading 
the collapse to OSD black holes is a non-generic subset of space of {M(r), b(r)}. The shortcoming of the tangential pressure 
case being that is not wholly physically satisfactory. Therefore to investigate the perfect fluid case, we moved to a 'bigger' space 
X x X x X, since the initial data set comprises of {M(r, Vi),p(r, Vi), b(r)}. Thus, in this space, the initial data set Iosd or 
Itbl or the union of both the sets, will become non-generic. Mathematically speaking, this set is meager or nowhere dense in 
the space X x X x X. This is proved by the study performed here. Thus, the initial configurations for the end states in the case 
of LTB or OSD models lie on the critical surface separating the two possible outcomes of collapse as discussed above. 

This analysis in fact shows how the structure of Einstein equations is very rich and complex, and how the introduction of 
pressures in the collapsing cloud opens up a lot of new possibilities that, while showing many interesting dynamical behaviours, 
do not rule out either of the two possible final outcomes. 

There are physical reasons, however, for the perfect fluid model to be subdued to the choice of an equation of state and there 
is also increasing evidence that such an equation of state cannot hold during the whole duration of the dynamical collapse. In 
fact there are indications that as the collapsing matter approaches the singularity large negative pressures arise, thus making the 
equation of state relating density to pressures depart from usual well-known equations of stellar equilibrium. Nevertheless the 
study of similar scenarios with linear or polytropic equations of state can give insights in the initial stages of collapse of a star. 
All this is very important from astrophysical point of view where still little is known of the processes that happen towards the 
very end of the life of a star, when in a catastrophic supernova explosion the outer layers are expelled and the inner core collapses 
under its own gravity. 

As we mentioned, due to the intrinsic complexity of Einstein equations for perfect fluid collapse, it is generally possible 
to solve the system of equations only under some simplifying assumptions (like the choice of a specific mass function), and 
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only close to the center of the cloud. The indications provided by the present analysis are then a first step towards a better 
understanding of what happens in the last stages of the complete gravitational collapse of a realistic massive body. 

Furthermore the above formalism could possibly be used as the framework upon which to develop possible numerical sim- 
ulations of gravitational collapse. As seen in the comoving frame, the positivity of \i (or X2)) is the necessary and sufficient 
condition for the singularity to be visible, at least locally. Numerical models of a collapsing star made of a perfect fluid with 
a polytropic equation of state (or a varying equation of state that takes into account the phase transitions that occur in matter 
under strong gravitational fields), with the addition of rotation and possibly electromagnetic field might help us better understand 
whether the inner ultradense region that forms at the center of the collapsing cloud when the apparent horizon is delayed, might 
be visible globally and have some effects on the outside universe. Many numerical models that describe dynamical evolutions 
leading to the formation of black holes exist both in gravitational collapse as in merger of compact objects such as neutron stars 
(see e.g. B73I1 . J7J]). but a fully comprehensive picture of what happens in the final moments of the life of a star is still far away. 

Close to the formation of the singularity, gravitationally repulsive effects, possibly due to some quantum gravitational correc- 
tions, are likely to take place. If such phenomena can interact with the outer layers of the collapsing cloud they might create a 
window to the physics of high gravitational fields whose effects might be visible to faraway observers. This scenario might in 
turn imply the visibility of the Planck scale physics or new physics close to the singularity, the presence of a quantum wall that 
might cause shock- waves from within the Schwarzschild radius that might give rise to different type of emissions and explosions, 
with photons or high energy particles escaping from the ultradense region. Collisions of particles with arbitrarily high center of 
mass energy near the Cauchy horizon might also happen [ 75I1 . 

Overall, the analysis led over the past few years seems to suggest that the Oppenheimer-Snyder-Datt scenario is indeed too 
restrictive to account for the richness of realistic dynamical models in general relativity. The occurrence of naked singularities in 
gravitational collapse appears to be a well established fact and a lot of intriguing new physics might arise from the future study 
of more detailed collapse models. 
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